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Abstract 

Acousto-Optical Coherence Tomography (AOCT) is variant of Acousto Optic Imaging (called also ultrasonic 
modulation imaging) that makes possible to get z resolution with acoustic and optic Continuous Wave 
(CW) beams. We describe here theoretically the AOCT effect, and we show that the Acousto Optic "tagged 
photons" remains coherent if they are generated within a specific z region of the sample. We quantify the z 
selectivity for both the "tagged photon" field, and for the M. Lesaffre et al. photorefractive signal. 

OCIS codes : 170.3660, 110.7050, 110.7170, 160.5320, 170.3880 



1 Introduction 

Acousto-optic imaging (AOI) [1, 2, 3] is a technique 
that couples ultrasounds and light in order to reveal 
the local optical contrast of absorbing and/or scatter- 
ing objects embedded within thick and highly scat- 
tering media, like human breast tissues. 

First experiments used fast single detectors to 
record the modulation of the optical signal at the 
US frequency [4, 5, 6, 1, 7]. But, since the phase of 
the modulation is different for each grain of speckle, 
the detector can only process one grain of speckle. To 
increase the optical etendue of detection, Leveque et 
al. [8] have developed a camera detection technique 



that processes many speckles in parallel. This tech- 
nique has been pulled to the photon shot noise limit 
by Gross et al. [9] using a holographic heterodyne 
technique [10] able to detect photons with optimal 
sensitivity [11, 12]. Since the US attenuation is low in 
tissues, the tagged photons are generated along the 
US propagation z axis with a nearly constant rate. 
This means that in a continuous regime of the US, 
the AO techniques give nearly no information on the 
location of the embedded objects along the z axis. 
To get such z information, Wang et al. [13], have de- 
veloped a US frequency chirp technique with a single 
detector, which has been extended to camera detec- 
tion [14, 15]. Unfortunately, these chirp techniques 
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cannot be used in living tissues, because the phase 
of light decorrelates very fast in them, since half fre- 
quency linewidth of light that travels through 4 cm 
of living breast tissue is about 1.5 kHz [16]. This 
phase decorrelation drastically lowers the detection 
efficiency, since the detection bandwidth is approxi- 
mately equal to the camera frame rate. The band- 
width is then much narrower than the width of the 
scattered photon frequency spectrum, and most of 
the tagged photons are undetected. It is still possible 
to increase the detection bandwidth by using a faster 
camera, but in such systems this generally means that 
a smaller number of pixels should be used, and the 
optical etendue of detection decreases accordingly. 
To perform selective detection of the tagged photons 
with high optical etendue, narrow band incoherent 
detection techniques have been proposed. For exam- 
ple, Li et al. select the tagged photon by spectral 
holeburning [17, 18], while Rousseau uses a confocal 
Fabry-Perot interferometer [19]. This last experiment 
[19] benefits from a powerful long pulse laser, whose 
duration (0.5 ms) matches the 1.5 kHz signal band- 
width. Another way to get a detection bandwidth 
comparable with the signal bandwidth while keeping 
a large optical etendue, detection schemes involving 
photorefractive (PR) crystals have been proposed. 
Murray et al. [20, 21] use a PR crystal sensitive at 
532 nm to select the untagged photons, which are 
detected by a single avalanche photodiode . In this 
case, the weight of the tagged photon signal is mea- 
sured indirectly by using the conservation law of the 
total number of photons (tagged + untagged) [22]. 
Ramaz et al. selectively detect either the tagged or 
the untagged photons [23]. The Ramaz technique is 
also able to measure in situ the photorefractive writ- 
ing time (tpr), which characterizes the detection fre- 
quency bandwidth [24]. 

In order to get information on the location of the 
object along the z axis, acoustic pulses can be used. 
The method has been extensively used both with sin- 
gle detectors [25, 26], cameras [3], PR crystals with- 
out [20, 21, 27, 28], or with long pulse laser [29]. 
Nevertheless, reaching a millimetric resolution with 
US pulses requires a typical duty cycle of 1%, corre- 
sponding to the exploration length within the sample 
(~ 10 cm) and the desired resolution (<~ 1 mm). This 



is problematic regarding the very small quantity of 
light that emerges from a clinical sample, since weak 
duty cycle yields low signal and poor Signal-to-Noise 
Ratio (SNR). When US pulses are used with photodi- 
ode detection, the SNR becomes lower, since fast pho- 
todetectors mean larger electronic noise. In a recent 
publication, Lesaffre et al. [30] overcome the duty 
cycle problem, and get z resolution with CW light 
and ultrasound by applying a random phase modula- 
tion on both the optical illumination and US beam. 
This so called Acousto-Optical Coherence Tomogra- 
phy (AOCT) technique is then demonstrated with 
photorefractive detection of the tagged photons. 

Whatever the method used to obtain an axial res- 
olution, the acousto-optic signal is sensitive to the 
quantity of photons tagged by the ultrasound : as 
shown in many previous studies, a strong absorber 
("zero" transmission") within the US field will induce 
an important drop on the signal [30, 1, 11, 18, 19, 22]. 
It has been shown more recently that a small "quasi- 
transparent" inclusion having a scattering coefficient 
(/j,' s = 10 cm -1 ) different from the host matrix (fi' s = 
7 cm -1 ) can give a contrast in the acousto-optic sig- 
nal [31]. In both cases, and to our knowledge, no 
quantitative measurements of this contrast have been 
performed as a function of the absorption coefficient 
nor the transport mean free path length I* . 

In the present paper, we will describe the AOCT 
effect theoretically. We show that the tagged pho- 
tons remain coherent if they are generated within a 
specific z region of the sample. We will quantify the 
z selectivity for both the tagged photon field, and for 
the tagged photon photorefractive signal as detected 
by Lesaffre et al. [30]. The theoretical results we get 
here will be compared with experiment in another 
publication. 

2 Theory of the Acousto Op- 
tic Coherent Tomography 
(AOCT). 

The theoretical description of the Acousto Optic Co- 
herent Tomography cannot be simply extrapolated 
from the theory made previously [22] to describe the 
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path ( „ r = (x,y,:) 

Figure 1: Light scattering along the travel path I, 
which involves scattering event located in r/, m , where 
m is the scattering event index. 

photorefractive detection of the UltraSound Modu- 
lated photons (USM) . Since we make tomography, we 
cannot consider that the USM photons are globally 
generated by the modulation of the length of a travel 
path. We must make a finer analysis by describing 
how the USM photons are locally generated within 
each specific region of the sample. 

2.1 The generation of the "tagged 
photons" 

Let us call Ep and E$ the fields coming into and out 
of the sample. Consider the point (x, y) located after 
the sample output interface. E$ is a quasi monochro- 
matic wave at the frequency ojq of the incoming laser. 
Let's introduce the complex field amplitude E_ P and 
E_ s defined as: 

E P (t) = ${E P e juJot } (1) 

E s (t) = ft{E s e^ ot } (2) 

where 5ft is the real part operator. Eg results from the 
sum (or the interference) of the field components E§j 
scattered through the sample along many travel paths 
I from input plane (z — 0) to the detector. Moreover, 
as illustrated by Fig.l, each travel path I can be de- 
composed in a succession of scattering events (Z,m) 



located in r^ m where m is the scattering events index. 

E s (t) = Y,E s ,i(t) (3) 
i 

= 5ft j^] a £ P e f-o(t-» l /c)| 

where / is the travel path index, and s; the corre- 
sponding effective travel path length. The length 
si is the product of the travel path length by the 
medium refractive index n. To simplify the discussion 
we consider that the field amplitude a is the same for 
all the travel paths. All travel paths have the same 
weight aE P , but different field phases: e~ J "° Si / c = 
e -j 2Trs,/A Since the travel path lengthes are large 
with respect to the optical wavelength A, the factor 
e -j27rsi/A j g ranc i onl Summing over the travel paths, 
one gets a speckle outgoing field. 

2.1.1 The ultrasonic field of pressure 

Let us now apply a CW (Continuous Wave) ultra- 
sonic (US) wave to the system by using an ultrasonic 
piezoelectric (PZT) device. The PZT transducer ex- 
citation voltage is: 

UpzT{t)=^{UpzT^ ust } (4) 

where U_pzt l& the complex amplitude of Upzt- Like 
in experiments, we consider here linear conditions 
where the acoustic pressure Pus is proportional to 
the excitation voltage. By this way, we get in any 
point r of the sample: 

P us (r, t) = A(r) U PZT (t - z/c us ) (5) 

where cus 1S the sound velocity in the sample, and 
z I cus the time delay from the US emission point (the 
PZT) to the zone of coordinate z that is considered. 
Let us introduce the US pressure complex amplitude 

£-US : 

Pus(r,t) - 5ft {P us (r) e^} (6) 

with 

P us (r) = A(r) Upz T e~^'^ (7) 

The pressure P_us IS periodic with respect to the 
US propagation axis, the period being Xus — 
2-kcus/uus- 
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2.1.2 The acousto optic modulation 

Because of the US beam the scatterers vibrate. More- 
over, the sample refractive index is modulated. These 
two effects yield a modulation of the length s; of the 
travel paths of the photons that are scattered by the 
medium (where I is the travel path index) at the US 
frequency ujjjs'- 



8l (t) =s lfi -$l{5s l e^ us '} 



(8) 



where <5s z is the complex amplitude of the modulation 
of the travel path I. We get from Eq.3: 



Es (t) 



??E a ^-p e 



j wo(t-si,o/c) 



(9) 



x exp 



Jwus t 



Let us introduce the complex amplitude Ss l m of the 
m th scatterer contribution to the travel path modu- 
lation, whose modulus and phase are /3/ jm and (f>i tTn 
respectively. 



(10) 



The sample outgoing field E s (t) is then modulated 
by the US at frequency lj us . 



Es(t) 



E P e 



j ui (t-si, /c) 



(11) 



x exp 



We get in Eq.ll: 



exp 



(13) 



The field Es{t) diffused by the sample becomes: 



aJ^Ep e*"o(t-«i.o/c) 



(14) 



The field Es(t) diffused by the sample is the sum of 
a main component i?s )U)0 (t), whose frequency is uj , 
with the two sideband components i?g )U , ±1 (t), whose 
frequencies are lu±± = loo ± ojus- 

E s (t) = E StU0 (t) + E StU1 (t) + Es^At) (15) 

Let us introduce ^5. Wo and E_ Smj±i , which are slow 
varying with time. 



= K {E s<Uo eM3"ot)} 
E s , u±1 (t) = » {E StU±i exp(jLJ ±1 t)} 
We get from Eq.14: 



(16) 



E StU0 (t) =»| aJ^Epe 



j u (t-si, /c) 



(17) 



E s ^{t) + Es^-At) = ® { aE P e j " ot 



(18) 



2.1.3 The tagging of the scattered photons 

In typical experiments, the vibration amplitude is 



E 

l.m 



.2tt/3 ; , 



+ C.C. 



much lower than the optical wavelength A = 2ttc/u> : where c c means the comp lex conjugate. We thus 

(19) 



for example, the vibration amplitude is 60 nm for 1 have for E 



MPa acoustic pressure at ujus — 2 MHz. We can 

then make the hypothesis of a weak acousto optic E Su±i (t) — aE P ^ 
modulation: ; 



^Esfe^^i « i 



(12) 



j e -j27rsi,o/A x 



E 



4 



Here, the main component Es jUo does not depend on 
the travel path modulation (Eq.17), while the mod- 
ulated components Es,u ±1 do. Moreover, whatever 
the modulation mechanism is: displacement of the 
scatterers or modulation of the refractive index, /3; jm 
is directly related to the acoustic pressure Pus (rj, m ) 
at the scatterer location r;. m . 

Note that the phases <pi^ m and (pi. m > of two scatter- 
ing events m and m' of the same path I are partially 
correlated according to the position of the associated 
diffusers r;. m and r;. m ', and according to the physical 
effect at the origin of the modulation. 

For the displacement of the scatterers, the phases 
4>i t m is related to the projection q z of the scattering 
wave vector q ; m along the US propagation direction 
(i.e z) with q, m = kj m - k Lm (where k Lm and k' l m 
are the wave vectors of the photon before and af- 
ter the scattering event l,m). The phases <fii, m and 
4>i_ m ' are not correlated, since q z may change of sign 
from one scattering event (/, m) to the next (I, m + 1) 
within the same path I. 

For the modulation of the refractive index, <f>i^ m is 
mainly related to the US phase. In a typical experi- 
ment the scattering length l s is about 0.1 mm, while 
the US wavelength Ays is about 1 mm (0.75 mm for 
ujjjs — 2 MHz). This means that (j)^ m and 4>i,m' are 
correlated, if the scattering events (I, m) and (l,m') 
are close together (\m — m'\ <a few units), and mi- 
correlated if not. 

This partial coherence allows us to use the acousto- 
optical modulation in scattering media. However, all 
the scatterers r; m in the acoustic column contribute 
to the tagged photons field £j (f), Thus on the 
acoustic column, the information is not localized. So 
it is necessary to use a complementary technique in 
order to obtain an axial z resolution. 

2.2 The axial resolution along z 

To obtain an axial resolution along z, Lesaffre et al. 
[30] have used Acousto Optic Coherent Tomography 
(AOCT). This technique is based on the control of 
the acoustic and optical coherence lengths using a 
random phase modulation on the acoustic and optical 
arms. 



2.2.1 The AOCT random modulation of the 
optical and acoustical field phases. 
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Figure 2: Order of magnitude of the various times, 
l/wo: optical period; 1/ujjjs- acoustic period; T<t>: 
acousto-optical correlation time; r: time averag- 
ing characteristic time; T mo d'- characteristic time of 
the modulation H(t); rp^: photorefractive time; r c : 
Lock-in integration time. 

The incoming optical field and the PZT excitation 
voltage are now: 

E' P {t) = {Epe^+^W >} (20) 

U' pzT (t) = Sft \U_ PZT e 3(^st-^us(t) ) j ( 21 ) 

where ipp et ipus are random phase modulations ap- 
plied to the optical incoming beam E' p and to the 
PZT that generates the US beam. Since we consider 
the effect of a random phase modulation, fields are 
noted Ep, E' s . The random phases ipp(t) an d tpus{t) 
are supposed to be fully correlated as follow: 

ipp (t) = ipus (t - za/c us ) (22) 

where zq/cus — is a fixed temporal delay which 
determines the z selected zone z ~ z$. 

To simplify the discussion we will consider here, 
like in [30], that the US phase ipus is randomly drawn 
every X$ to be or 7r with equal probability. The 
optical phase ipp follows the same random phase law 
than ipus, but the phase is dealyed in time by Zq/cus- 
The incoming complex field is: 

& P {t)=Epe^ p ^ (23) 



and the US excitation (U_'pzt)> an d US pressure 
(Pjjs) complex amplitudes are: 



U_' PZT (t) 
E! US {r,t) 



u 



PZT 



jipus(t) 



(24) 



A( r ) Hpzt e 



ju us zt/c us e jipus(t-z/cus) 



Thus we eliminate the fast varying components of 
E^g (t) w hich anyway will have no effect on the 
final signal. To be complete let's define here the tem- 
poral average operator ( ) T : 



2.2.2 The "tagged photons" field. 

By making the calculations leading to Eq.19 with the 
random phases ipp and ipus, we get the tagged pho- 
tons complex amplitude E! S u : 



T J t ' = 



t'=t+r/2 



t>=t-r/2 



(....)df 



(28) 



E's i0J±1 (t) = aE P J2 
i 



-j27rs !i0 /A 



(25) 



±3<t 



±A,m(t) 



A 



where the phase Vi,m ; which depends on time t, and 
on location z t _ m of the m t/l scatterer along the axis 
z, is defined by: 

1pl,m(t) = 1pp(t) + TpUS {t ~ Zl,m/c US ) ■ (26) 

Because of the random phase jumps, which occur ev- 
ery T^, the complex field E^ S u]±i (t) varies with a char- 
acteristic time T^, while, in absence of random mod- 
ulation, the field E_ s ^ ±i does not depend on time. 
In the following, we will detect the field £j (t) by 
photorefractive effect on a crystal. 

We must notice that all the photorefractive detec- 
tion E' Suj±i (t) processes occur on times much larger 
than T#. 

• The writing of photorefractive signal on the crys- 
tal occurs in a time rpp 3> T^. 

• To get a modulated signal for the Lock-In am- 
plifier, the phase of the US will be modulated at 
a frequency oj mod = l/T mod with T mod > Tq. 

• The extraction of the modulated signal modu- 
lated at LOmod will be made via Lock-In with an 
integration time r c ^> T$. 

So one can replace in the following the field E[ S u±i (t) 
by its temporal average (E' 1 (t)) r over the charac- 
teristic time r chosen such as (see Fig. 2): 



Ta, < T -C Tp_R, T mod , T c 



(27) 



The temporal average of the tagged photon field over 
the characteristic time r is then: 



(E' s , u±1 (t))r = aE P 



J e 



-j2irsi, /\ 



('Z—lLHL, 



A 



x (e 



±A,m(t) 



(29) 



As we can see on Eq.29, ipi im acts on the temporal av- 
erage (E' StU±i (t)) T only through (e ± ^ i .">W) r , which 
depends only on the location along z of the scatterer 
of indexes l,m, i.e. on z^ m . From Eq.22 and Eq.26, 
we have 



*,m(t) - 

(e ± ^.-(*)) T ~ 1 



(30) 



ipl,m(i) — 0,7r randomly 
( e =y*.™(*)) T ~ (31) 

for the scatterer i,m located in (Eq.30) and out 
(Eq.31) the selected zone z ~ zo respectively. 

To characterize this 2; selection mechanism in a 
more quantitative way, let us define the time correla- 
tion function: 



£i(*) = 



(gp(t)E£(t + 9)) T 

(\E' P (t)\\ 



(32) 



= l e j^us{t) e -j^us{t+e)^ 



{\Ulus{t)\\ 

In the case of 0,7r random phase jumps considered 
here, g (6) is a triangular function that corresponds 



6 




Figure 3: Plot of correlation function g {&) (a), and 
its square |.g 1 (6')| 2 (b). The horizontal axis Units is 
cither 0/7$ (for time correlation), or (z — zq)/Az with 
Az — cjjs 7$ (f° r z resolution). 
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Figure 4: Principle of the photorefractive detection. 



to the convolution of two rectangles of width Tpy^. 
The correlation function g^ff) is plotted on Fig. 3 (a). 

The field (E_ Sul ) can be expressed as a function 
of g{. 



& Wi (<)>t = €pE 

><E 



J e 



-j'27rs ii0 /A 



(33) 



CUS 



A 



Let us note here that the second member of Eq. 33 
does not depend on time. It means that, when we 
apply the random modulations of phase, the field 
B j (t) reaches, after a brief transitory regime, a 
stationary regime in which the slow field components 
do not depend on time any more. 

Furthermore, the results of the calculations do not 
depend on r as soon as the condition T^, C t < 
(tpr, T mo d, t c ) of Eq.27 is fulfilled. So one should 
write: {g s>u±1 (t)) r = i^S^J- 

2.3 The photorefractive detection of 
the tagged photons 

We will now consider the photorefractive detection 
of the tagged photons in order to quantify the z se- 
lection process for the photorefractive detected sig- 
nal (and not just for the tagged field E[ s u itself). 
The calculations we will make are similar to the ones 



made by Gross at al. [22], but in a slightly different 
context. 

2.3.1 The detection principle 

The principle of the photorefractive detection is il- 
lustrated on Fig. 4. The signal Eg, the wave front of 
which is distorted, is collected in a photorefractive 
crystal. A reference beam Ep, considered as plane 
wave, which is also called pump beam, interferes with 
it within the crystal. By photorefractive effect, the 
interferogram grooves a hologram corresponding to a 
weak modulation Sn(r) of the local refractive index 
within the volume of the crystal. This effect having a 
finite response time Tpp, only the static component 
of the interferogram contributes to the recording of 
the hologram. 

To simplify the analysis, we will consider the de- 
tection of tagged photons of the +1 sideband at 
uj\ = ujq + lojjs- So we will shift the beam refer- 
ence frequency by ojjjs — uji —ujq in order to perform 
the photorefractive detection at frequency oji . Let us 
introduce the complex amplitude E_ R of the refer- 
ence field. 



E R {t) = 3? [E Rtl 



(34) 



The photorefractive effect selects, in the signal field 

^5 = E S^ + E 's^i + E S,u-i> the field com P° Iient 
E' s . The reference beam is then diffracted by the 
holographic grating grooved within the crystal yield- 
ing the field E' D , whose wavefront is the same for 
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E'sui- At ^ ne ex ^ °f the crystal, one gets then 
both the transmitted signal beam E' s , and the beam 
diffracted by the crystal, which will be noted E' D . 



2.3.2 The reference field diffracted by the 
crystal E' D 

Let us introduce the complex amplitude of the 
diffracted field defined by 



transient components. By making the transforma- 
tion t' — > (i — t'), in the integral J dt' of Eq.3, we 
obtain: 



— P - a (y-yi)/2 



Els, ui {y,t) 

Jo 



(38) 
f) G(y,t') 



(35) 



Let us call y — y 1 and y = y 2 the crystal en- 
trance and exit planes respectively, and t = the 
origin of time, when no photorefractive hologram is 
recorded. Let us consider that the reference field 
E_ R is constant. Within the crystal, the signal field 
Els oji (j/i < V < V2, t) can be written as a function of 
the entrance field E StU)1 {yi,t) [32, 33] 



We can notice that the hologram is written with de- 
layed time t — t' with a delay t' varying from zero 
to some tpji. Also let us note that the first term of 
Eq. 38 corresponds to the signal wavefront that is 
transmitted by the crystal. Let us call EL its field 
and E[ T its complex amplitude (the index T means 
here transmitted). The second term corresponds to 
the reference field that is diffracted by the crystal we 
will note E' D (where the index D means diffracted), 
and E[ D for the complex amplitude. We can write: 



E 



<*(v-Vi)/2 



W,u»i(».*) = 

(»!.*) + 



dt'El s ^{ yi ,t')G{y,t 



(36 fsM) 

f T (y,t) 



E' T (y,t)+E' D (y,t) (39) 
E' T ^ (V, *) + E' TiU0 (y, t) + E' TtU+i {y, t) 

E' D , u+1 (y,t) 



where, under conditions of weak recording efficiency 
and weak absorption, the transfer function G(y, t) 
can be written as [34] 



G(y,t) 



i(y-yi) 



e T ra 



TPR 



(37) 



Here, t pr is the photorefractive response time and 7 
the photorefractive gain. Equation 36 is established 
for two plane waves, but it can be generalized to dis- 
torted wavefront by decomposing the wavefront in 
plane waves. Several approximations are made to 
establish this equation: (i) the reference beam is a 
monochromatic wave with constant frequency uj\ (i.e. 
it is not temporarily modulated), (ii) it is not per- 
turbed by the recording of the hologram although it 
can be attenuated by the crystal, and (iii) its power 
is larger than the signal beam one. 

In the AOCT experiment [30], the tagged photon 
signal is modulated at a frequency w mod of some kHz, 
and Lock-in detection is performed. So we are inter- 
ested in the low-frequency evolution of E_ s . So 
we can replace the t' = lower limit of the inte- 
gral / dt' of Eq. 36 by t' = —00, by neglecting the 



gs^AV't) = E! T ^{y,t)+El DtUJi {y,t) (40) 

&r, u Av,t) = e-^y-y^ 2 E' s ^( yi ,t) 

x / dt'E! s ^{y u t-t')G(y,t') 
Jo 

Note that since the photorefractive effect selects the 
field components of frequency LO\, the diffracted field 
E' D exhibit in Eq.39 a single frequency component 

E' T 



D. 



OJ+l 



2.4 The acousto optic signal detected 
by a large area photodiode 

We consider that the signal is detected by a photo- 
diode of large area located near the cristal exit plane 
y = j/2- The photodiode signal Spd is equal to the 
integral of \El s \ 2 over its area. We get from Eq.39. 

S PD {t) = c.c.+ (41) 

J dx J dz \E' s (x,y 2 ,z,t)\ 2 
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S PD (t) =c.c.+ 
J dx j dz | \E? T (x, 2/2, z, t)\ 2 

+ \E' D {x,y 2 ,z,t)\ 

+ {E! T (x,y 2 ,z,t) E'£(x,y 2 ,z,t))} 



Because the acousto-optical interaction does not 
modify the total number of photons, i.e. tagged 
+ untagged photons, the term |-E T | 2 in Eq.42 does 
not depend on the acoustic modulation. Further- 
more, because the gain is supposed to be low, i.e. 
Vpr = 7(2/2 — Vi) < lj the term \E^ D \ 2 can be ne- 
glected in front of the crossed term E[ T E[ D . There- 
fore we can only consider the crossed term, which is 
the product of the diffracted field E' D , which builts up 
with the characteristic time rpp, and the transmit- 
ted field E' T , which can vary quickly. The photodiode 
modulated signal S PD is thus: 



S' PD {t) = c.c. + (43) 
J dx J dz E! TiUi {x,y 2 ,z,t)E , £ UJi (x,y 2 ,z,t) 

From Eq.40, we get 

S' PD (t) = c.c.+ (44) 
e -«(»-f0 j dx J dzE' StUi (x,y u z,t) 

x / dt' El ui (x, yi ,z,t-t') G*(y 2 ,t') 
Jo 



To keep a certain universality, we write G* although 
G is supposed to be real. We can then develop E' s Wi 
by summing up all the paths (index I) and scatter- 
ing events (index m) contributions by using Eq.33. 



(42) Averaging over a time r, we get: 



(S' PD (t)) T =c.c. + e- a ^- v ^\aE p \ 2 J dx j dz 
^2j e -j2™,, /\ x ^ 2lrf3l ' m n ( z ° Z ZLm ~\ 



J^dt' G*(y 2 ,t')x (^2-je 2 i™".° 



A - 1 V vus 



\ " 27T^;>,m' * ( ZQ - Zy 

^-f A - 1 I v us 



(45) 



The equation 45 illustrates the complexity of the cal- 
culation of the signal. It involves a double summation 
over the optical paths ( i.e. ,,), a double summa- 
tion over the scattering events (i.e. J2 m m'), a spatial 
integral over the photodiode area (i.e. J J dx dz), and 
a temporal integral over the delay t' (i.e. J dt'). 

To simplify this equation, let us consider first the 
integral over the photodiode area J J dx dz. Every 
point (x,z) of the photodiode selects paths I and I', 
which finishes in (x, z). For the corresponding paths, 
the phase factor e"^ 271 " Si / A is totally random from a 
route to the next one. So we can limit the summation 
over I and V to the terms I — I'. The equation 45 
becomes then: 



(S' PD (t)) T = ex. + e 



-a(V2-yi) 



2naE p 



EE 



l m 



ZQ - Zl >r 
CUS 



Zq - Zl t 
cus 



J dx J dz 

X) 

dt'G*(y 2 ,t') 

J{4>l,m-4>l, m ') 



(46) 



To simplify this equation further, it is necessary to 
study the mutual coherence of the phases <pi !m and 
4>i_ m i that corresponds to two different scattering 
events m and m' of the same path I. According to 
the position of the scatterers, and according to the 
acousto-optical modulation mechanism, these phases 
are correlated or not. Nevertheless, when the two 
events (m and m!) occur in two z coordinates z^ m 
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and zi !m > separated by more than an acoustic wave- 
length X us (i.e. z (jTO / —zi tm \ > Xus), the phases </>; iTO 
and <pi t m' are weakly correlated. We can then write: 



to X$: the scattering events that contribute to the 
signal must be within the g (z) selected region. 



(S' PD (t)h 



-a(y 2 -yi) 



= c.c. + 
2naE r 



X 



EE E 

( m TO ' with \zi, m - 



cj(4>l,-n 



(47) 

CO 

J dt'G*(y 2 ,t') J dx j dz 

ZQ — Zl^ m 



i |<Aus 



ZQ - Zl, 
CUS 



9, 



cus 



Since the magnitude of the acoustic pressure vary 
weakly over Xus, we have (3^ m ~ /3{ >m / for z; iTO - 
^,to' < A(/s- Moreover, the random modulation 
of phases is chosen in such a way that the char- 
acteristic length cus is larger than the acoustic 
wavelength Xjjs (i- e - T<f, cus > ^us- This implies 
that g l ((z - zi m )/c us ) ~ g l ((z - z^ m ,)/c us ) for 
zi to — to' < Ays. Therefore we obtain: 



+e 



-oi{v2-yi 



2waE r , 



X 



(48) 

dt'G*{y 2 ,t') j dx j dz 



EEa 



I m 

X 



cus 



e^' 

to' with |2; )Tn — 2^ m / |<At/s 



Here, the term 



selects the zone 



g 1 (zo - z Lm /c us ) 

of imaging. 

The z resolution one can expect is roughly equal 
to the half-width of g^z), i.e. to 0.5 x CusT® (see 
Fig. 3). The expected resolution is thus about 7.5 mm 
for T* = 10 ^s (20 US periods at uj us = 2 MHz), and 
1.1 mm for 7$ = 1.5 fis (3 US periods at u us = 2 
MHz). The AOCT published experimental results 
[30] correspond to T$ ~ 2 /j,s. To improve the z reso- 
lution, one must thus decrease T$. The acousto optic 
signal decreases accordingly, since it is proportional 



2.5 The Lock-in detection of the 
acousto optical signal 

2.5.1 The modulation of the signal at uj mo d 

Note that (SpD(t)) T given by Eq.48 is invariant with 
time. So the tagged photons photorefractive signal 
Spd is a CW component, which adds to the total flow 
of transmitted light. To detect Spd more efficiently 
with a Lock-in, AOCT adds an extra modulation of 
the ultrasonic wave. Like in the AOCT experiment 
[30] , we will consider here an asymmetric to tt phase 
modulation Hjj g (t) at frequency Lo mo d = 27r /T, mo( i ~ 
3 kHz, with duty cycle < r < 1: 



Hus{t) = +1 pour 
Hus(t) = -1 pour 



< t/T mod < <49) 
r < t/T mod < 1 



We consider that the modulation period T mo d is 
very large compared to the correlation time T^, 
but smaller than the photorefractive time rpp, i.e. 
(2tt/lous) < T<j, < Tmod < t pr . In practice, we typi- 
cally use T moc i ~ 100/zs (see Fig. 2). The modulation 
is applied according 

U' PZT (t) -> U! PZT (t) = H us (t) U! PZT {t) (50) 

The US signal is denoted C/'pzt an d the fields are 
denoted E' p , Eg and so on. The complex amplitude 
of the tagged photon field is now 

E!s, Ul (*) -»• E! S , U1 (*) = H(t) E! StUJi (t) (51) 

where H(t) = Hus{t—zv/cus)- Similarly with Eq. 40, 
the diffracted complex amplitude becomes 

/>oo 

El D ^{y,t) = e- a ^' 2 / dt'[H{t-t') 

Jo 

xE!s,«M,t-t') G(yX)} (52) 
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2.5.2 The modulated acousto-optical signal the signal on a large area photodiode thus becomes. 



The signal from the photodiode given by equation 44 
becomes 



(S PD (t)) r ~c.c + 



-a(y2-Vi) 



2naE p 



S PD (t) = c.c. 
+ e~^ jdxjdz 

/ dt' H(t-t)E'*s, Ul < 
Jo 



(53) xfdxfdzJ2Y.fr 



(56) 

H(t)(l-2r) 7(2/2-2/1) 
2 

I Zo — Zl.m \ 

£l 



H(t)E' s ^(x,y 2 ,z,t) 
{x,y 2 ,z,t-t') G*(y 2 ,t') 



I m 
m' /\z l-m -z l ,\<\ us 



CUS 



"By this way, the photodiode signal is modulated fol- 
lowing H{t). 

By making the calculation leading to Eq. 48 from to Equation 56, which does not depend on t can be 
Eq. 44 with the additional modulation H us , we get. slightly simplified as following: 



(S PD (t)) T = C.C.+ 



(54) 



(S PD (t)) = (1 - 2r)H(t)\E p \ Z j(y 2 - yi )e 



-a{v2-yi) 



-ct(y2-Vi) 



2-KClE r 



H{t) jdt'H(t-t')G*(y 2 ,t') J dX J dZ ^^2 



Pi, m 



9i 



zp — z^ m 
Cus 



JdxJ dzJ2Y, 



I m 



Pi 



zo — z^ m 
cus 



(e j (^-^')+c.c) 

m' with \zi rrn — z l m /\<Xus 



E 



' l\r-l, m -z lm i\<\us 



(57) 



3 Conclusion 



Since we have consider T mod < t pr , the integration The mains results of the P a P er are Ec l- 33 and E( l- 57 - 
over t' can be simplified, and we obtain from Eq.37. Equation 33 shows that the tagged photon field 

< s s, w± i>. which is 



/ dt'H(t-t')G*(y 2 ,t') = 
Jo 



Tpr 



Tpr 



dt'H(t - t') 



X 


[/ 




.JO 




(1- 



dt'G*(y 2 ,t') 



• proportional to the amplitude of the optical field 

(55) 

E P injected in the scattering medium, 

• proportional to the acoustic power delivered by 
the PZT via the term fif >m 



g^zo - z Um /c us ). 

This means that for a modulation H{t) faster than The random phase modulation creates along the 
t P r, the photorefractive recorded hologram is pro- acoustic column a zone of coherence located near 
portional to the average {H{t)) TpR . The asymmetric z ~ zq = cjjsO. The tagged photon signal from 
nature of the modulation H (t) yields non-zero pho- this zone adds up coherently, and can be further 
torefractive grating. The modulated component of detected. For the tagged photon field, the random 
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phase jump z selection is quantified by the factor 
g (zo — zi, m /cjjs)i which is the correlation product 
of a rectangle of width T$. This correlation product, 
which has triangular shape, is plotted on Fig. 3(a) as 
a function of time, in X$ Units, or as a function of 
the scatterer z relative coordinate (i.e. z — zi, m ), in 
Az = cjjsTq, Units. 

On the other hand, Eq.57 shows that the acousto- 
optical modulated signal on the large area photodiode 
surface (SpD(t)) r is 

• proportional to the optical intensity \Ep\ 2 in- 
jected in the scattering medium, 

• proportional to the surface of the photodiode, 
i.e. J dx J dz, 

• proportional to the acoustic power delivered by 
the PZT via the term /3f m 

• proportional to (1 — 2r)H(t) where H(t) is the 
additional time modulation, whose duty cycle is 
r. Because of this well controlled time modu- 
lation, the photodiode signal can be Lock-in de- 
tected at the frequency oj m od with an integration 
time t c > T mod = 2ir/uj mod . 

• and proportional to the square of the correlation 

2 

function g v i.e. to g^zo - z^ m /c us ) 

For the photodiode signal, the random phase 
jump z selection is quantified by the factor 

2 

g x (zo — z^ m /cus) , which is the square of the cor- 
relation product g^ This IgJ 2 factor is plotted on 
Fig. 3(b). One must note also that the summation 

over m! of the phases factors e J ( < ^' m ~ < ^> m ') , which can 
be limited to \z^ m — zi_ m > \ < Xus, describes here the 
effect of partial coherence of the successive scattering 
events within a given travel path I. The correspond- 
ing proportionality factor does not depend on the US 
modulation, and does not provide any z selection. 

In this paper, we have described the AOCT ef- 
fect theoretically. We show that the tagged photons 
remain coherent if they are generated within a se- 
lected z region of the sample, and we have quanti- 
fied this z selection effect for both the tagged pho- 
ton field E' s , and the photorefractive photodiode 



signal S PD (t). These theoretical results will be com- 
pared with experiment in another publication. 
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